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the  energetics  of  photoionization  are  calculated  for  transparent  and 
absorbing  solvents:  electronic  polarization,  London  dispersion  and  Born 
repulsion  energies  for  a  discrete  model  of  coordinated  solvent  molecules 
in  the  inner-sphere  solvation  shell  of  anions  and  cations;  electronic 
polarization  of  the  outer-sphere  region  for  a  continuous  medium  model. 
The  losses  resulting  from  the  rapid  variation  of  the  ionic  field  for  an 
absorbing  solvent  are  calculated  for  the  inner-  and  outer-sphere  regions, 
respectively,  from  a  discrete  model  and  a  continuous  medium.  Damping  of 
the  ionic  field  resulting  from  solvent  absorption  is  negligible.  The 
theory  is  applied  to  aqueous  solutions  in  the  7  to  10.4  eV  range  of 
photon  energies  by  using  dielectric  data  from  reflectance  spectroscopy  of 
liquid  water.  Experimental  dispersion  spectra  for  photoelectron  emission 
have  the  shape  predicted  by  theory  and  display  all  the  extrema  at  the 
photon  energies  of  the  calculated  curves.  The  very  pronounced  effect  of 
ionic  strength  on  the  balance  between  inner-  and  outer-sphere 
contributions  predicted  by  theory  (inner-outer  sphere  splitting)  is  fully 
confirmed  by  experiment.  Dispersion  spectra  of  inorganic  ions  in  the 
range  of  each  of  the  two  absorption  bands  of  liquid  water  (maxima  at  ca. 
8.2  and  10.0  eV)  therefore  exhibit  a  double  maximum  for  normal  dispersion 
and  a  double  minimum  for  anomalous  dispersion  (12  extrema  between  7.2  and 
10.4  eV).  Specific  effect  of  the  nature  of  anions  is  evident  above  9.0 
eV  in  inner-outer  sphere  splitting.  The  present  study  provides  a  way  of 
probing  the  response  of  liquids  and  solutions  to  the  rapidly  varying 
intense  ionic  field  resulting  from  the  process  of  photoionization. 
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Nonequilibrium  electronic  polarization  of  the  solvent  in  photoionization 
Paul  Delahay  and  Andrew  Dziedzic 

Department  of  Chemistry,  New  York  University,  New  York,  New  York  10003 


(Received  ) 

The  energetics  of  photoionization  in  condensed  phases  includes  a 
significant  contribution  from  nonequilibrium  processes  arising  from  dielectric 
dispersion  of  the  solvent  at  the  prevailing  photon  energy.  The  solvent  is 
polarized  by  the  varying  electric  field  caused  by  the  change  of  ionic  valence 
as  a  result  of  photoionization.  This  ionic  field  varies  in  a  time  interval 
determined  by  the  frequency  of  incident  radiation.  The  following 
contributions  from  nonequilibrium  processes  to  the  energetics  of 
photoionization  are  calculated  for  transparent  and  absorbing  solvents: 
electronic  polarization,  London  dispersion  and  Born  repulsion  energies  for  a 
discrete  model  of  coordinated  solvent  molecules  in  the  inner-sphere  solvation 
shell  of  anions  and  cations;  electronic  polarization  of  the  outer-sphere 
region  for  a  continuous  medium  model.  The  losses  resulting  from  the  rapid 
variation  of  the  ionic  field  for  an  absorbing  solvent  are  calculated  for  the 
inner-  and  outer-sphere  regions,  respectively,  from  a  discrete  model  and  a 
continuous  medium.  Damping  of  the  ionic  field  resulting  from  solvent 
absorption  is  negligible.  The  theory  is  applied  to  aqueous  solutions  in  the  7 
to  10.4  eV  range  of  photon  energies  by  using  dielectric  data  from  reflectance 
spectroscopy  of  liquid  water.  Experimental  dispersion  spectra  for 
photoelectron  emission  have  the  shape  predicted  by  theory  and  display  all  the 
extrema  at  the  photon  energies  of  the  calculated  curves.  The  very  pronounced 
effect  of  ionic  strength  on  the  balance  between  inner-  ind  outer-sphere 
contributions  predicted  by  theory  (inner-outer  sphere  splitting)  is  fully 
confirmed  by  experiment.  Dispersion  spectra  of  inorganic  ions  in  the  range  of 


each  of  the  two  absorption  bands  of  liquid  water  (maxima  at  ca.  8.2  and  10.0 
eV)  therefore  exhibit  a  double  maximum  for  normal  dispersion  and  a  double 
minimum  for  anomalous  dispersion  (12  extrema  between  7.2  and  10.4  eV). 

Specific  effect  of  the  nature  of  anions  is  evident  above  9.0  eV  in  inner-outer 
sphere  splitting.  The  present  study  provides  a  way  of  probing  the  response  of 
liquids  and  solutions  to  the  rapidly  varying  intense  ionic  field  resulting 
from  the  process  of  photoionization. 

I.  INTRODUCTION 

Dielectric  dispersion  of  the  solvent  was  shown  recently*  to  affect  the 

energetics  of  photoionization  of  ions  and  molecules  in  solution.  This  effect 
1  2 

was  interpreted  ’  by  assuming  that  the  energetics  of  solvation  of  the 
species  being  photoionized  is  determined  by  the  solvent  polarizability  a  and 
optical  dielectric  constant  eop  at  the  photon  energy  E  at  which 
photoionization  occurs  rather  than  by  the  limiting  value  a0  and  e°  . 
Alternatively,  one  may  assume  that  the  change  of  ionic  valence  of  the 
substance  being  photoionized  occurs  on  a  time  scale  determined  by  the 
frequency  of  the  incident  radiation.  The  response  of  the  medium  to  the 
rapidly  varying  intense  ionic  field  is  determined  under  these  conditions  by  a 
and  at  the  photon  energy  E  and  not  by  the  limiting  values  a0  and 
e°p*  The  dispersion  effect  according  to  this  model  results  from 
nonequilibrium  processes  involving  electronic  polarization  of  the  solvent, 
London  dispersion  and  Born  repulsion.  Energy  losses  resulting  from  the  rapid 
variations  of  the  ionic  field  must  also  be  considered  for  absorbing  'solvents. 
This  model  is  treated  quantitatively  and  tested  experimentally  in  the  present 
paper.  Differences  in  the  solvent  configuration  around  anions  and  cations  are 
taken  into  account  following  a  comment  by  Hush^  on  the  earlier  treatment. 


The  selection  of  the  proper  dielectric  constant  of  the  solvent  in  the 
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present  treatment  follows  from  the  three-step  model  used  in  the 
interpretation  of  photoelectron  emission  by  liquids  *  :  (i)  photoionization, 

(ii)  transfer  of  quasifree  electrons  to  the  liquid-gas  interface,  and  (iii) 
interfacial  barrier  problem.  Dispersion  affects  only  photoionization,  and  the 
dielectric  constant  eQp  is  selected  accordingly,  as  noted  above.  The  proper 
dielectric  constant  (or  range  of  dielectric  constants}  in  the  ion-electron 
separation  of  the  above  step  (ii)  is  determined  by  the  kinetic  energy  of  the 
quasifree  electrons  generated  by  photoionization.  The  present  paper  is  not 
concerned  with  this  problem  since  ion-electron  separation  and  transfer  are 
embodied  in  the  emission  law  accounted,  for  instance,  by  the  theory  of  Refs.  5 
and  6. 

II.  NONEQUILIBRIUM  PROCESSES:  TRANSPARENT  SOLVENTS 
A.  Theoretical  model 

7  8 

The  fairly  standard  treatment  of  ionic  solvation  ’  will  be  adapted  to 

the  present  problem  with  the  additional  consideration  of  the  different  solvent 

configurations  around  inorganic  anions  and  cations.  Electrically  neutral 

species  can  be  treated  in  terms  of  the  continuous  medium  model  of  Sec.  IIC 

without  any  particular  difficulty.  The  volume  around  the  ion  being 

photoionized  is  divided  into  two  regions,  (i)  The  inner-sphere  solvation 

shell  consists  of  N  solvent  molecules  coordinated  to  the  ion.  The  solvent 

molecules  are  treated  as  point  dipoles  at  a  distance  r  =  r  +  r  from 

r  r  o  c  w 

the  ionic  charge,  where  r  and  r  are  the  crystallographic  radii  of  the 

c  w 

ion  and  solvent,  respectively,  (ii)  The  outer-sphere  region  outside  the 
sphere  of  radius  a  =  r  +  2r  is  treated  as  a  continuous  medium. 


The  nonequi librium  contribution  from  electronic  polarization  in  the 
present  treatment  is  determined  by  the  change  az  =  1  in  the  ionic  valence  and 
thus  is  independent  of  the  value  of  z  as  in  the  case  of  nuclear  reorganization 
of  the  outer-sphere  region  (Marcus,^  Hush*®).  A  multipole  expansion** 
of  the  potential  of  the  ion  will  be  used  as  in  the  calculation  of  energies  of 
ionic  solvation.  The  polarizability  of  the  solvent  is  assumed  to  be 
isotropic.  This  is  nearly  the  case  for  water  considered  below  in  the 
application  of  the  theory.  Calculated  energies  differ  by  only  ca.  1  percent 
if  anisotropy  is  taken  into  account  on  the  basis  of  the  polarizability  data  of 
Ref.  12.  The  induced  dipole  vector  is  collinear  with  the  electric  field 
vector  in  the  central  field  of  the  ionic  charge  if  one  assumes  isotropic 
polarizability  of  the  solvent.  Nonlinear  effects  are  not  considered. 


B.  Inner-sphere  solvation  shell 

Only  the  terms  depending  on  a  and  eop  of  the  solvent  and  the  terms  for 

London  dispersion  and  Born  repulsion  are  retained  in  the  expression  for  the 

7  8 

free  energy  of  ionic  solvation.  ’  We  calculate  first  the  change  of  energy 
aPin  for  the  inner-sphere  solvation  shell  for  az  *  1  and  the  change  of 
polarizability  from  a  to  a0.  One  has 

Api-  =  aUIND  +  AU(Pae)  +  AU(Pap)  +  AU(PaPa)  +  AU(Paq) 


in 


+  aUl(S,s)  +  aUl(s,s)  +  All 


REP’ 


(1) 


where  UjNq  represents  the  work  required  for  the  formation  of  induced  dipoles 
p^  in  the  inner-sphere  solvation  shell;  the  next  four  terms  are  the  energies 
of  interaction  of  induced  dipoles  with  the  ionic  charge  (pae),  dipoles 
(p  p),  induced  dipoles  (p  p  ),  quadrupoles  (p  q);  U.  (S.s)  and 

a  a  a  a  L 

UL(s,s)  are  the  energies  for  solute-solvent  and  solvent-solvent  London 

dispersion,  respectively;  Up£p  corresponds  to  solute-solvent  and 

solvent-solvent  Born  repulsion. 


Explicit  forms  are: 


AUIN0  =  (N/2)[(p°)2/a°  -  p2/a] 

(2) 

Al)(pae)  =  -  (Ne/r2)(p°  -  pj 

(3) 

AU(pop)  =  2A(p/r2)(p°  -  pjcos  e 

(4) 

AU(PaPa)  -  (A/r2)[(p°)2  -  p2] 

(5) 

AU(poq)  =  -  B{0/rJ)(p°  -  pa) 

(6) 

aUl(S,s)  =  -  (3N/2r®)[IsI/(Is  +  I)]as(a°  -  a) 

(7) 

AUL(S,S)  =  -  (CI/r06)[(a°)2  -  a2] 

(8) 

aUREP  =  ■  (1/x)t2AU(Pae)  +  3[AU(pap)  +  AU(poPa)] 

+  4AU(pQq)  +  6[aUl(S,s)  +  aUl(s,s)]} 

(9) 

where  N  is  the  number  of  solvent  molecules  coordinated  to  the  species  being 
photoionized;  e,  the  absolute  value  of  the  electronic  charge;  p,  the  permanent 
dipole  moment  of  the  solvent;  q,  the  quadrupole  moment  of  the  solvent 
(Appendix);  a,  the  angle  between  the  induced  and  permanent  dipoles  of  the 
solvent  in  the  inner-sphere  solvation  shell;  I$  and  I,  the  gas-phase 
ionization  energies  of  the  solute  and  solvent,  respectively;  a$,  the 
polarizability  of  the  solute  taken  to  be  independent  of  photon  energy  in  the 
absence  of  dispersion  data.  The  coefficients  A,  B  and  C  were  calculated  by 
Buckingham^  by  considering  solvent  pairwise  interactions:  A  =  2.296  and 
7.114;  B  =  1.722  and  5.336;  C  *  0.2373  and  1.160  for  tetrahedral  and 


octahedral  coordination,  respectively.  The  Born  repulsion  energy  UREp 

Q 

obtained  according  to  Morf  and  Simon  by  minimizing  with  respect  to  rQ 


is 

the 


total  contribution  from  the  inner-sphere  shell  to  the  solvation  energy.  The 


exponent  x  will  be  taken  to  be  12  and  8,  respectively,  for  inorganic 


The  induced  dipole  moment  p°  is  obtained  according  to  Morf  and 
8  o 

Simon  by  minimizing  with  respect  to  p  the  total  energy  depending  on 

a 

p°.  Thus, 

p°  *  a°(Ner^  -  2AprQcos  6  +  B0)/rQ(Nr^  +  2Aa°)  (10) 

The  same  equation  with  o  instead  of  a0  is  applied  to  p  .  Equation  (10) 

a 

reduces  to  the  familiar  equation  p°  =  a°F,  where  F  =  e/rjj  is  the 

14 

field,  if  the  pairwise  interactions  are  neglected.  The  sum  of  the 
energies  in  Eqs.  (2)  and  (3)  in  that  case  is  one-half  of  the  value  on  the 
r.h.s.  of  Eq.  (3). 

The  evaluation  of  the  angle  8  in  Eq.  (4)  and  0  (Appendix  A)  requires 

information  about  the  structure  of  the  solvent  in  the  inner-sphere  shell. 

Only  inorganic  anions  and  cations  in  aqueous  solution  will  be  considered 

here.  The  anion,  hydrogen  and  oxygen  are  collinear  in  the  simplest 

case*5-^7  and  the  water  dipole  vector  makes  an  angle  8  =  52.23°  with  the 

electric  field  vector  of  the  ionic  charge.  The  water  molecules  are  assumed  to 

13 

rotate  freely  about  the  anion-hydrogen-oxygen  axis.  Conversely,  the 
cation  and  water  molecule  are  assumed  to  be  in  the  same  plane  in  the  simplest 
case  with  the  cation-oxygen  axis  coinciding  with  the  electric  field  vector  of 
the  ionic  charge^-^7  (8  =  0).  The  water  molecule  is  assumed  to  rotate 
freely  about  the  cation-oxygen  axis.7  Tilting  of  the  plane  of  the  water 
molecule  with  respect  to  the  cation-oxygen  axis,  which  was  reported^-^7  for 
high  electrolyte  concentrations  (>  1.5  M) ,  will  not  be  considered  here  but 
could  easily  be  introduced  in  the  treatment.  To  summarize,  the  limiting  values 
8  =  52.23°  and  8=0  will  be  used  for  anions  and  cations,  respectively,  in  the 
application  of  the  theory  to  aqueous  solutions. 


7 


C.  Screening  by  the  ionic  atmosphere 

The  change  in  the  free  energy  for  electronic  polarization  of  the 
outer-sphere  region  is  treated  for  a  continuous-medium  model.  One  has 

AP0ut  -  (e^aJn.oV1  -  <U> 

in  the  absence  of  screening  of  the  ions  by  the  ionic  atmosphere.  A  correction 


must  be  made  for  the  screening.  The  free  energy  contribution  AGatm  from  the 
ionic  atmosphere  to  the  free  energy  of  electronic  polarization  of  the 


outer-sphere  region  is, 

aGatm  <e2'2‘op>K,(1  +  «>•  <12> 

to  the  approximation  of  the  Debye-Huckel  theory.  There  is  the  optical 

dielectric  constant,  a  is  the  radius  of  the  inner-sphere  taken  to  be  equal  to 

r„  +  r, .  (Sec.  IIA),  and  <  is  the  Debye-Hlkkel  reciprocal  length  measured 

L  W  1  1  1 

from  a.  The  quantity  <  is  calculated  for  the  static  dielectric  constant  es 
of  the  solvent  because  the  ionic  atmosphere  is  "frozen"  during  the  rapid 
variations  of  the  ionic  field  resulting  from  photoionization.  Equation  (12) 

■JO 

is  taken  directly  from  the  Debye-Huckel  theory  except  that  eQp  is 

written  in  the  denominator  instead  of  £$.  A  correction  identical  to  Eq. 

19 

(12)  applies  to  nuclear  reorganization  in  optical  electron  transfer  except 

that  e _  appears  in  the  denominator  instead  of  e  . 
a  op 

It  is  seen  from  Eq.  (12)  that  the  correction  for  the  ionic  atmosphere  is 
obtained  by  writing  (e  /2R) (1/ eQp)  with  a  radius  R  =  a  +  1 l<  instead  of 
R  =  a.  The  energy  SaPQut  therefore  is  obtained  instead  of  APQut  of  Eq.  (11), 
the  screening  factor  being, 

S  =  [a-1  -  (a  +  l/<)'1]/a“1 

=  1/(1  +  <a).  (13) 

Thus,  the  value  S  =  1  corresponds  to  the  limiting  case  of  no  screening 
approximated  for  very  dilute  solutions.  Conversely,  the  value  S  =  0  pertains 


to  total  screening  of  the  outer-sphere  response.  One  has,  for  instance,  for  a 
0.1  M  aqueous  solution  of  a  1-1  electrolyte,  K  ~  6.S  A  denoting  the 
thickness  of  the  ionic  atmosphere  outside  the  inner-sphere  of  radius  a  =  4.76 
A;  consequently,  S  -  0.6.  Thus,  the  effect  of  screening  by  the  ionic 
atmosphere  in  aqueous  solutions  can  be  very  significant  for  electronic 
polarization  whereas  it  is  minor  for  nuclear  reorganization  in  optical 
electron  transfer  (a  few  hundredths  of  electronvolt  at  most  for  a  concentrated 
solution) . 

Expressions  for  the  ionic  screening  of  the  discrete  inner-sphere  are 
complicated  but  may  be  closely  approximated  by  a  simple  extension  of  the 
continuum  outer-sphere  treatment.  For  solutions  of  high  ionic  strength, 
counterions  may  penetrate  the  inner-sphere;  thus,  the  Debye-Huckel  length 
<-1,  measured  from  the  radius  a,  is  negative  denoting  this  penetration.  The 
screening  factor  S  is  then  negative  allowing  for  the  partial  cancellation  of 
the  inner-sphere  response  by  that  calculated  for  the  screened  outer-sphere 
(see  Sec.  IV).  For  example,  in  a  2  M  aqueous  solution  of  a  2-1  electrolyte, 
the  mean  separation  of  ionic  centers  is  ~  6.5  A.  Given  an  effective  common 

1  O 

ionic  radius  r£  =  2  A,  the  distance  from  one  of  the  ions  to  the  center 
of  the  other  is  4.5  A.  Thus,  the  inner-sphere  is  penetrated  to  a  depth  of 
4.76  -  4.50  =  0.26  A.  as  defined  above  is  then  ^-0.3  A  and 
consequently  S  is  ca.  -0.1. 

Formula  (13)  is  intended  as  a  guide  in  the  interpretation  of  screening 
rather  than  a  quantitative  equation.  Application  of  the  Debye-Huckel  theory 
outside  the  range  of  sufficiently  dilute  solutions  is  tentative,  but  the 
original  model  can  be  improved  (e.g.,  the  recent  treatment  in  Ref.  20).  It  is 
noted,  in  this  respect,  that  departure  from  the  theory  does  not  arise  because 
<  *  in  concentrated  solutions  is  shorter  than  ionic  radii  since  the  ionic 


atmosphere  is  always  outside  the  sphere  of  of  closest  approach  for  the  ions 
the  radius  of  which  is  approximately  the  sum  of  the  anionic  and  cationic 
radi i . 


III.  NONEQUILIBRIUM  PROCESSES:  ABSORBING  SOLVENTS 

A.  Real  part  of  energies  for  inner-  and  outer-sphere  regions 

The  dielectric  constant  of  the  solvent  in  the  range  of  absorption  is  a 
complex  quantity  -  ie-,*  The  energy  density  within  the  dielectric  is 
ejF2/8ir,  where  F  =  (e/r2)(e2  +  e2)”^2  is  the  magnitude  of  the 

2 

ionic  field.  Integrating  the  energy  density  for  the  volume  element  dV  =  4irr  dr 
between  r  =  a  and  r  =  one  obtains  the  outer-sphere  contribution 

*Pout  -  (e2/2a)C(c°p)-1  -  .j/t.2  *  •!)].  (14) 

Thus,  Eq.  (14)  is  simply  the  real  part  of  aPout  of  Eq.  (11)  written  for  a 
complex  dielectric  constant. 

Similarly,  calculation  of  aP^  requires  the  value  of  the  real  part 

of  the  complex  polarizabil ity  -ic^.  The  values  of  and  a?  are 

obtained,  at  least  approximately,  from  the  Lorenz-Lorentz  equation  for  an 
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absorbing  solvent 

(ej  -  1  -  ifigJ/Uj  +  2  -  ie£)  =  c(°i  “  i°2)»  (15) 

where  ?  =  4ttNa6/3M,  being  the  Avogadro  number,  s  the  density  and  M  the 

molecular  weight  of  the  solvent.  Equation  (15)  yields 

C ( ~  l)(ej  +  2)  +  e|]/[(El  +  2)2  +  e|]  =  (16) 

e2/[(e1  +  2)2  +  e2]  =  (C/3)a2  (17) 

The  energy  AP-n  is  computed  from  Eqs.  (1)  to  (9)  by  using  from  Eq.  (16) 

instead  of  a.  Equations  (2)  to  (6)  are  written  in  complex  form  in  terms  of 

o,  -  ia0  and  p  -  ip  (Sec.  1 1 1 B ) .  The  sum  of  these  complex 

energies  is  minimized  with  respect  to  P„  _  ip  ,  and  the  real  and 

1  a2 


imaginary  parts  of  induced  dipole  are  obtained  accordingly.  Thus, 

Pa  -  +  <*2)  ^1^2°1^^2  +  K^a^)  +  (Kgag)  3 

PQ  *  ®2^l^2^^2  +  ^3al  )  +  (^3a2^  ^ 

where 


(18) 

(19) 


K1  =  Nero  =  2Apro  cos  6  +  B 
K2  -  Nr« 

k3  -  2flro 


(20) 

(21) 

(22) 
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Use  of  the  corrected  Lorenz-Lorentz  equation  yielded  values  of 
and  (*2  slightly  higher  than  those  given  by  Eqs.  (16)  and  (17)  (cf.  also  Sec. 
V). 


B.  Imaginary  part  of  energies  for  inner-  and  outer-sphere  regions 

The  imaginary  part  e2  of  the  complex  dielectric  constant  can  be 
interpreted  classically  in  terms  of  the  optical  conductivity  a  =  iDe2/4ir  at 
the  photon  energy  *fl<u .  This  interpretation,  holds  even  if  the  dielectric  does 
not  exhibit  any  conductivity  from  free  charges,  and  the  imaginary  part  e? 
arises  solely  from  bound-bound  transitions.  (The  contribution  from  ionic 
conductivity  for  the  systems  studied  in  the  present  work  is  negligible.)  The 
energy  loss  corresponding  to  absorption  will  be  calculated  for  the  outer-  and 
inner-sphere  regions.  The  loss  for  the  outer-sphere  region  will  be  calculated 
first  since  a  continuous  medium  treatment  is  directly  applicable. 

The  conduction  current  density  for  the  field  F  is  aF.  The  corresponding 
energy  loss  in  the  dielectric  is  determined  by  the  variation  of  F  during  the 
time  interval  for  which  the  loss  is  calculated.  This  time  interval  is  taken 
to  be  equal  to  the  time  uncertainty  corresponding  to  the  photon  energy 


E  =  Tiu)  at  which  photoionization  occurs.  (See  formulation  of  the  time-energy 
uncertainty  in  Ref.  22.)  The  ionic  valence  cannot  be  specified  during  the 


interval  o>  and  consequently  the  loss  is  calculated  by  equating  F  to  the 
change  in  ionic  field  from  an  ionic  charge  z  to  z  +  1  as  a  result  of 

photoionization.  The  loss  per  unit  of  time  and  unit  volume  therefore  is 

?  21  2 
cF  .  This  result  differs  from  the  classical  value  oF  /2  obtained  from 

j  aEdE.  The  experimental  results  on  the  effect  of  ionic  atmosphere 

J°  2  2 
in  Sec.  VI  strongly  support  the  use  of  aF  instead  of  cF  /2. 

-1  -1  2 

The  energy  loss  per  unit  volume  in  the  time  u>  is  w  oF  or 

(e2/4n)F  .  The  ionic  field  in  the  absence  of  screening  by  the  ionic 
22  2  -1/2 

atmosphere  is  (e/r  )(e^  +  e£)  ,  and  the  volume  element  dV  is 

2 

4wr  dr.  The  energy  loss  in  the  outer-sphere  region  is 
Lout  ■  '2  / 

a. 

=  (e2/a)e2/(e2  +  e2).  (23) 

Equation  (23)  is  identical,  except  for  the  factor  1/2,  to  the  negative 
imaginary  part  of  the  complex  polarization  which  is  associated  with  the  real 
part  of  the  polarization  given  by  Eq.  (14).  The  correction  for  the  ionic 
atmosphere  is  obtained  just  as  in  Sec.  IIC  by  introducing  the  radius  a  +  1/* 
instead  of  a  in  Eq.  (23).  The  screening  factor  S  of  Eq.  (13)  therefore  is 
also  applicable  to  L  t. 

The  loss  in  the  inner-sphere  L.p  is  calculated  analogously  to  LQut  by 
summing  the  negative  imaginary  parts  of  the  complex  components  of  APir). 

This  is  justified  by  noting  that  the  conduction  current  density  per  molecule 
within  a  time  oT*  is  the  displacement  of  charge  out-of-phase  with  the  real 
part  of  the  ionic  field  due  to  optical  conductivity.  This  out-of-phase  charge 
displacement  is  alternatively  given  by  the  formally  defined  dipole  p 

a2  * 

The  dipole  p  defined  in  this  way  interacts  with  the  permanent  and 
“2 

induced  moments  of  the  multipole  expansion  to  give  inner-sphere  energy 
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losses.  In  this  way,  the  inner-sphere  interaction  energies  are  damped  by  the 
optical  conduction  of  the  medium.  Thus, 

Lin  *  L(0IND^  +  L(e<^  +  L^p°^  +  L^paa^  +  +  LREP  (24) 

where  *- ( 0 j ND )  rePresents  t*ie  damping  of  the  induction  energy  over  the 

interval  w-*  arising  from  the  optical  conductivity  of  the  solvent. 

Likewise,  L(ea),  L(pa),  L(pqo),  L(q,a)  represent  the  decrease  of  the 

interaction  energies  (or  losses)  for  the  charge,  dipole,  induced  dipole  and 

quadrupole  terms  of  Eqs.  (2)  to  (6),  respectively.  The  quantity  L(oJND)  is 

the  imaginary  contribution  from  the  Born  repulsion.  This  term  is  included 

because  the  minimization  of  the  interaction  energies  with  respect  to  nuclear 

coordinates,  which  leads  to  Eq.  (9),  minimizes  conductive  losses  as  well. 

London  dispersion  is  not  applicable  to  the  loss  calculation,  however.  One  has 


(cf.  Sec.  IIB) , 

L(W  ■  -  P l2)  -  ♦  “s’2  <26> 

L(eo)  .  (N/2)(e/rP)p  (26) 

o  a2 

L(po)  =  -  2A(pp  lr3) cos  a  (27) 

02  u 

L(p  o)  =  -  2Ap  p  /r3  (28) 

Qt  0^2  U 

L(qcr)  =  B0p  /rj  (29) 

u 

LREp  -  (l/x){2L(eo)  +  3[L(pa)  +  L(poo)]  +  4L(qo)}  (30) 


where  p  and  the  quantities  A,  B,  p,  e,  x  are  defined  in  Sec.  IIB. 
a2 

IV.  TOTAL  ENERGY  FOR  NONEQUILIBRIUM  PROCESSES 

The  total  contribution  W  from  nonequilibrium  processes  to  the  exclusion  of 
nuclear  reorganization  will  be  obtained.  The  real  and  imaginary  components  of 
W  are 


P  ■  lPin  *  lPout 


L  "  Lin  +  Louf 


(31) 

(32) 
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Thus, 

W  =  (P2  +  l2)112  (33) 

The  variations  of  APin,  AP0Ut>  Lin»  Lout  and  W  with  Photon  energy 
for  S  =  1  are  displayed  in  Fig.  1  for  liquid  water  for  octahedral  coordination 
of  the  solvent  around  an  anion.  The  shapes  of  the  curves  of  Fig.  1  are 
determined  by  dielectric  dispersion  and  the  absorption  bands  of  liquid  water 
at  ca.  8.2  and  10.0  eV.  This  point  is  amplified  in  Sec.  V  and  VI.  Results 
for  tetrahedral  coordination  around  an  anion  are  the  same  except  that  aPin 
makes  a  contribution  to  W  which  is  ca.  18  percent  lower  than  in  Fig.  1. 

Results  for  cations  are  similar  to  those  for  anions  for  the  same  coordination. 

Data  on  and  used  in  Fig.  1  were  those  recently  recalculated  by 

Painter  from  earlier  results  on  reflectance  spectroscopy  of  liquid 

24 

water  except  for  The  earlier  values  of  in  the  6.8  to  8.8  eV 
interval  were  used  because  the  cut-off  =  0  below  7.6  eV  in  Refs.  25  and 
26  does  not  agree  with  the  data  from  other  workers  summarized  in  Ref.  27.  The 
difference  is  really  minor  in  the  calculation  of  L  but  not  in  that  of  the 
derivative  dL/dE  needed  for  comparison  with  experiment.  The  quadrupole  moment 
components  used  in  the  preparation  of  Fig.  1  were  taken  from  Ref.  28. 

The  effect  of  screening  by  the  ionic  atmosphere  on  the  variations  of  the 
energy  W  of  Eq.  (33)  with  the  photon  energy  E  is  shown  in  Fig.  2  for  S  =  0.3 
to  S  =  -0.1.  The  shape  of  the  W  against  E  curve  hardly  changes  between  S  =  1 
and  0.3,  but  W  at  a  given  E  decreases  with  increasing  screening.  Figure  2,  in 
fact,  shows  that  screening  can  drastically  reduce  the  magnitude  of  W.  Whereas 
the  outer-sphere  contributions  APQut  and  LQut  become  dominant  for  S  £,0.2, 

one  has  only  N  .  <iP?n  ♦  for  the  case  of  s  .  0.  Negative 

S  values  even  further  reduce  W  since  the  inner-sphere  response  is  decreased  by 

counterion  penetration  (Sec.  IIC). 
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V.  CALCULATED  DISPERSION  SPECTRA 

1  2 

The  present  theory  was  tested  experimentally  by  using  data  *  on 

photoelectron  emission  by  aqueous  solutions  of  inorganic  anions  and  cations  in 

the  7  to  10.4  eV  range  of  photon  energy.  Verification  of  the  theory  requires 

plots  of  -dW/dE  against  the  photon  energy  E.  This  will  be  shown  first.  The 

emission  yield  Y  in  emission  experiments  is  defined  as  the  number  of  collected 

photoelectrons  per  incident  photon.  One  has  Y  =  K(E  -  Et)p,  where  E  is 

the  photon  energy,  Et  the  threshold  energy,  K  a  quantity  independent  of  E 

and  p  =  2,  5/2  or  3  depending^  on  the  range  of  E.  A  quadratic  dependence 
1  2 

was  found  suitable  ’  and,  in  any  event,  conclusions  to  be  drawn  below  are 
essentially  unaffected  by  the  choice  of  exponent. ^  One  has  in  view  of  the 
form  of  the  emission  law  for  Y, 

dY1/2/dE  =  K(1  -  dEt/dE)  (34) 

29  30 

The  threshold  energy  Et  is  given  by  * 

Et  -  aGcyde  *  R  *  “•  <35> 

where  AGCyC-|e  can  be  written  in  explicit  form  for  a  thermodynamic  cycle 

involving  the  initial  and  final  states  of  the  emission  process  (including  the 

surface  potential  contribution),  and  R  is  the  free  energy  for  nuclear 

reorganization.  The  quantities  AGc^c-|e  and  R  are  not  affected  by  dielectric 

dispersion  and  are  independent  of  E.  One  obtains  from  Eqs.  (34)  and  (35) 

dY1/2/dE  =  K[1  -  dW/dE].  (36) 

It  follows  from  Eq.  (36)  that  the  derivatives  dY^2/dE  and  -dW/dE  should 

have  the  same  functional  dependence  on  photon  energy.  The  curves  representing 

1/2 

the  variations  of  dY  / dE  and  -dW/dE  with  E  will  be  referred  to, 
respectively,  as  the  experimental  and  theoretical  dispersion  spectra 
From  Eq.  (33) 

dW/dE  =  [P/(P2  +  L2)1/2]dP/dE  +  [L/(P2  +  L2)1/2]dL/dE. 


(37) 


The  two  terms  on  the  r.h.s.  of  Eq.  (37)  were  computed  as  functions  of  photon 
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energy  by  using  binomial  smoothing  filters  to  obtain  the  derivatives  dP/dE 

and  dL/dE.  This  method  is  preferable  to  the  use  of  the  Savitsky-Golay 

1  32 

convolutes  previously  used  ’  to  calculate  the  derivatives.  It  should  be 

noted  in  this  respect  that  the. data  on  and  liquid  water  from  Ref.  23 

used  in  the  preparation  of  Fig.  1  were  available  at  0.2  eV  intervals  and  were 

interpolated  by  means  of  quadratic  and  quartic  polynomials.  Interpolation  was 

of  no  consequence  in  the  calculation  of  the  curves  of  Fig.  1  but  resulted  in 

the  elimination  or  at  least  damping  of  any  fine  structure  the  derivative 

curves  might  actually  have  in  the  0.2  eV  intervals.  Figure  3  shows  that  the 

terms  in  dP/dE  and  dL/dE  in  Eq.  (37)  make  comparable  contributions  to  -dW/dE 

but  that  their  extrema  occur  at  different  photon  energies. 

21 

The  use  of  the  corrected  Lorenz-Lorentz  equation  increases  the 
magnitude  of  J-dW/dt|  by  ca.  10  percent  but  does  not  change  the  functional 
dependence  of  -dW/dt  on  photon  energy. 

The  difference  in  the  orientation  of  the  water  molecules  in  the 
inner-sphere  shell  around  anions  and  cations  (Sec.  I I B)  was  investigated. 
Calculated  dispersion  spectra  for  cations  exhibit  practically  the  same 
functional  dependence  on  photon  energy  as  anions  but  the  amplitude  is 
increased,  e.g.,  by  ca.  50  percent  in  comparison  with  Fig.  3  for  the  following 
data:  N  =  6,  r  =  0.88  A,  e  =  0,  o  =  1.3  x  10'24  cm3,  I  =  26.4  eV,  x  =  8 
(data  corresponding  to  the  nonhydrated  V  in  Ref.  8). 

The  effect  of  screening  by  the  ionic  atmosphere  is  displayed  in  Fig.  4. 

The  shape  of  the  dispersion  spectrum  curve  is  not  sensitive  to  screening 
between  S  =  1  and  S  =  0.2,  but  a  major  change  in  shape  occurs  at  S  $  0.2. 

This  is  the  case  because  the  curves  corresponding  to  dP/dE  and  dL/dE  peak  at 
different  photon  energies  (Fig.  3).  The  curve  for  S  =  0  corresponds  to  the 
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case  in  which  the  dispersion  spectrum  is  solely  determined  by  the  inner-sphere 
contributions.  Conversely,  the  contributions  from  the  outer-sphere  region  are 
dominant  for  S  =  1.  The  photon  energies  at  the  extrema  of  the  dispersion 
spectra  calculated  for  S  =  1  and  S  =  0  are  listed  in  Table  I.  The  ranges  of  E 

<  7.75  eV  and  7.75  <  E  <  8.70  eV  correspond  respectively  to  normal  and 

anomalous  dispersion  for  the  first  absorption  band  of  liquid  water  (with 

maximum  at  ca.  8.2  eV).  A  maximum  between  the  two  minima  is  displayed  in  the 

normal  range,  and  conversely  a  minimum  between  two  maxima  is  observed  in  the 
anomalous  range.  This  pattern  of  extrema,  which  will  be  referred  to  as 
inner-outer  sphere  splitting,  is  repeated  in  the  normal  and  anomalous  ranges 
at  higher  photon  energies  for  the  second  absorption  band  of  liquid  water  (ca. 
10.0  eV). 

VI.  COMPARISON  OF  EXPERIMENTAL  AND  CALCULATED  DISPERSION  SPECTRA 

Experimental  and  calculated  dispersion  spectra  are  compared  in  Figs.  5  and 
6  for  1-1  (NaSCN)  and  1-2  ^828203)  electrolytes.  The  ordinate  of 
experimental  dispersion  spectra  is  obtained  in  arbitrary  units,  and 
consequently  the  curves  of  Figs.  5  and  6  were  normalized  at  the  maximum  (8.2 
to  8.6  eV)  and  minimum  (ca.  9.1  eV)  of  the  experimental  curves.  The 
calculated  curves  were  prepared  for  S  =  0.3  and  0.15  (Fig.  5)  and  S  =  0.0  and 
-0.1  (Fig.  6)  for  the  0.25  and  2  M  solutions,  respectively.  These  values  of  S 
are  significantly  different  at  high  ionic  strengths  from  the  ones  computed 
from  Eq.  (13)  with  measured  from  a:  0.41  (0.25  M)  and  -0.14  (2  M)  for  a 
1-1  electrolyte;  0.14  (0.25  M)  and  -0.40  (2  M)  for  a  1-2  electrolyte.  This  is 
not  surprising  since  an  accurate  calculation  of  from  the  Debye-Huckel 
theory  in  its  original  form,  as  applied  in  Eq.  (13),  is  not  possible  at  the 
concentrations  of  Figs.  5  and  6. 
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The  experimental  curves  for  the  2  M  solution  of  Figs.  5  and  6  have 

slightly  different  widths  than  the  calculated  dispersion  spectra.  Three 

possible  reasons  at  least  can  be  advanced  for  this  discrepancy:  the 

uncertainty  about  and  inherent  to  the  application  of  the 

33 

Kramers-Kronig  transformation  in  reflectance  spectroscopy  ;  the  uncertainty 
about  ^2  below  ca.  8.8  eV  (Sec.  IV);  the  approximation  of  the  theoretical 
model  in  which  spherical  symmetry  of  the  field  is  assumed  around  the  ion 
undergoing  photo  ionization. 

The  change  in  the  shape  of  dispersion  spectra  upon  increase  of  the  ionic 
strength  is  strikingly  apparent  in  Figs.  5  and  6.  Thus,  the  maximum  with 
dominant  outer-sphere  contribution  near  8.6  eV  is  prominent  for  the  0.25  M 
sodium  thiocyanate  solution  whereas  the  maximum  with  dominant  inner-sphere 
contribution  near  8.2  eV  is  considerably  enhanced  for  2  M  thiosulfate. 
Inner-outer  sphere  splitting  (cf.  Table  I)  in  the  8.0  to  8.8  eV  and  9.0  to  9.7 
eV  ranges  is  indeed  observed.  Figures  5  and  6  also  show  that  the  weak  maximum 
at  9.38  eV  in  the  calculated  curves  is  greatly  enhanced  in  the  experimental 
curves.  Inner-outer  sphere  splitting  in  the  9.0  to  9.7  eV  range  therefore  can 
be  strongly  influenced  by  the  nature  of  the  anion.  This  specific  effect  is 
displayed  by  the  17  anions  studied  in  Ref.  1.  Thus,  the  height  of  t(ie  9.37  eV 
peak  is  nearly  independent  of  concentration  (0.25  to  2  M)  for  some  anions 
(e.g.,  OH  ,  CO^-,  SCN  )  whereas  it  is  strongly  concentration-dependent 
for  others  (e.g.,  Br“,  I",  CIO^).  The  9.63  eV  minimum  is  also 
sensitive  to  the  nature  of  the  anions  as  is  evident  from  Figs.  5  and  6. 
Particularly  pronounced  minima  were  observed  in  Ref.  1  with  bromide  at  9.62  eV 
and  iodide  at  9.59  eV.  Greatly  enhanced  inner-outer  sphere  splitting  was 
recently  observed  for  phosphorus  hexafluoride  ion  (0.2  M  KPFg,  Et  ~  9.3  eV) 
with  the  normal  sequence  of  two  maxima  (10.02  and  10.35  eV)  and  a  minimum 


(10.18  eV).  This  enhancement  from  PFg  ion  completely  masked  the 
dispersion  spectrum  of  liquid  water  (E^  =  10.06  eV). 

The  photon  energies  at  the  extrema  of  the  experimental  dispersion  spectra 
of  the  17  anions  studied  in  Ref.  1  in  Table  I  compare  very  well  with  the 
calculated  values  for  -0.2,$  S  ,$  1.  It  is  concluded  from  this  agreement  and 
the  interpretation  in  Sec.  V  and  VI  that  the  theory  accounts  remarkably  well 
for  the  experimental  dispersion  spectra  and  the  pronounced  effect  of  ionic 
strength  on  their  shape. 

The  present  theory  is  supported  also  by  the  observed  shift  in  threshold 
energy  with  concentration  of  ionic  emitter  for  NaSCN,  NaN^,  Na^S^O^ 
and  Na2S20g  over  the  0.25  to  2.0  M  range.  Further  work  on  this  aspect 
will  be  the  subject  of  a  future  publication. 

VII.  CONCLUSIONS 

Photoionization  provides  a  way  of  changing  the  intense  ionic  field  in  a 
time  interval  corresponding  to  the  frequency  of  the  incident  radiation.  These 
rapid  variations  of  the  ionic  field  can  be  used  to  probe  the  response  of 
liquids  and  solutions.  Experimental  evidence  was  obtained  in  this  way  for 
nonequilibrium  electronic  processes  resulting  from  dielectric  dispersion  of 
the  solvent  in  photoionization.  Theory  predicts  and  experiment  confirms  that 
ionic  screening  can  alter  profoundly  the  balance  between  the  contributions 
from  inner-  and  outer-sphere  regions  to  the  energetics  of  the  nonequilibrium 
electronic  processes  (inner-outer  sphere  splitting). 

ACKNOWLEDGMENTS 

This  work  was  supported  by  the  Office  of  Naval  Research  and  the  National 
Science  Foundation.  We  thank  Professor  N.  S.  Hush  (Theoretical  Chemistry, 
Sydney)  for  pointing  out  to  us  the  need  to  consider  the  difference  in  solvent 
configuration  around  inorganic  anions  and  cations.  Professor  L.  R.  Painter 


19 


(Physics,  University  of  Tennessee)  kindly  communicated  her  unpublished 
recalculated  data  on  the  reflectance  spectroscopy  of  liquid  water. 


APPENDIX 

A.  Quadrupole  moment 

The  quadrupole  moment  of  Eq.  (6)  for  6  =  0  is' 


0=0  +0. 
a  b 


where 


0+0.  =20  -0  -  0  , 
a  b  zz  xx  yy’ 


the  z-axis  being  in  the  ion-oxygen  direction.  For  e  =  52.23°,  one  has 


9a  +  eb  *  2ezzsec  6  -  Vec  8  -  V 


Equation  (A3)  was  obtained  by  noting  that  the  water  molecule  rotates  about  the 
H-0  axis  for  e  =  52.23°.  The  quantities  and  ©b  therefore  were 
transformed  with  respect  to  the  new  axis  of  rotation. 


B.  Attenuation  of  the  field 

34 

The  electric  field  in  a  conducting,  isotropic  medium  is  attenuated  by 

the  factor  exp(-uikr/c) ,  where  k  is  the  absorption  coefficient,  r  the  depth  of 

penetration  and  c  the  velocity  of  light  in  vacuum.  The  absorption  coefficient 

2  2 

k  can  be  computed  from  =  n  -  k  and  =  2nk,  where  n  is  the  index 

9 

of  refraction.  The  field  in  the  present  case  is  proportional  to  r  without 

attenuation.  The  energy  without  attenuation  is  obtained  by  integrating 
-4  ? 

r  dV  from  a  to  <*>,  where  dV  =  4*r  dr.  With  attenuation,  one  integrates 
4irr  exp(-2wkr/c)dr  between  a  and  «>.  This  is  done  in  the  calculation  of 
APout  using  the  real  part  of  (e^  -  Uj)-*  as  noted  in  Sec.  III.  There 
is  no  attenuation  of  the  field  for  the  limiting  value  e°p.  One  obtains 
APout  =  (e2/2a)(c0°p)-1  -  (e2/2)El(^  +  e2)"1[a_1exp(-ca) 

+  Sin  sa  -  S2a  +  C3/2.2!  ...],  (A4) 

where  s  =  2<uk/c.  The  series  in  (A4)  converges  rapidly  in  the  present  case  and 


only  the  first  term  between  brackets  is  significant.  Equation  (A4)  reduces  to 
Eq.  (11)  for  5=0.  One  has^  for  liquid  water  at  8.2  eV,  k  =  0.19,  5  = 

•3  1 

7.5  x  10  cm  and  exp(-£a)  =  0.992.  Attenuation  of  the  field  thus 
changes  the  second  term  on  the  r.h.s.  of  Eq.  (A4)  by  only  *  1  percent. 
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TABLE  I.  Photon  energies  at  the  extrema  of  experimental  and  calculated 
dispersion  spectra  (S  =  0  and  1) 


Dispersion 


normal 


(<  7.75  eV) 


Extremum  Dominant 


inner 


mixed 


outer 


Calculated3  Experimental^ 


Contribution  Photon  Energy  Photon  Energy 


no  datum 


anomalous 


(7.75  to  8.70  eV) 


normal  min 

(8.70  to  9.70  eV)  max 


anomalous 
(>  9.7  eV) 


inner 


mixed 


outer 


inner 


mixet 


outer 


inner 


mixed 


outer 


10.11 

10.32 


8.16*0.06 

8.42*0.06 


8.65*0.03 


9.11*0.03 

9.37*0.06 

9.63*0.12 


9.86*0.04 


10.18 


10.35 


See  Fig.  4  for  7  4:  E  ^  10  eV. 

bData  for  0.5  M  VC  1 2  and  1  M  CrCl2  (7.30  and  7.41  eV)  and  0.05  M  K4Fe(CN)6 
(7.41  eV)  from  Ref.  2.  Average  values  with  standard  deviation  (8.16  ^  E  .<:  10  eV) 
for  1  M  solutions  of  17  inorganic  anions  from  Ref.  1.  Extrema  at  10.18  and 


10.35  eV  obtained  with  0.2  M  KPFg. 


Captions  to  Figures 


Fig.  1  Variation  of  energies  with  photon  energy  for  liquid  water  and  the 

solvent  configuration  around  a  typical  inorganic  anion:  aP^,  and 

a P  L  ,  for  the  inner-  and  outer-sphere  regions,  respectively,  W  the 

out  out 

total  energy  from  Eq.  (25).  Data:  N  =  6,  rc  =  2  A,  rw  =  1.38  A,  B  = 

52.23°,  p  =  1.855  debye,  a0  =  1.444  x  10-24  cm3,  a$  =  4  x  10~24 
cm3,  I  =  3  eV,  I  =  12.61  eV,  x  =  12.  Values  of  and  e2  from  Ref. 

23  (see  text).  No  screening  by  the  ionic  atmosphere  (S  =  1). 

Fig.  2.  Variations  of  the  total  energy  W  with  photon  energy  for  different 
values  of  the  screening  factor  S  (S  =  0  for  complete  outer-sphere  screening,  S 
=  -0.1  for  partial  screening  of  the  inner-sphere).  Same  data  as  in  Fig.  1. 

Fig.  3.  Variations  of  the  quantities  in  Eq.  (29)  with  photon  energy  for  the 
data  of  Fig.  1. 

Fig.  4.  Variations  of  -dW/dE  with  photon  energy  for  different  values  of  the 
screening  factor  S  (S  =  0  for  complete  screening  of  the  outer-sphere,  S  =  -0.1 
for  partial  screening  of  the  inner-sphere).  Same  data  as  in  Fig.  1. 

Ordinates  at  the  minimum  at  7.42  eV  from  S  =  -0.1  to  S  =  1.0:  -0.164,  -0.196, 

-0.229,  -0.263,  -0.297,  -0.541.  Ordinates  at  the  maximum  (photon  energy 
between  parentheses)  from  S  =  -0.1  to  S  =  1.0:  0.151  (8.20  eV),  0.151  (8.24 

eV),  0.162  (8.35  eV),  0.181  (8.41  eV),  0.216  (8.68  eV),  0.509  (8.68  eV). 

Fig.  5.  Comparison  of  experimental  dispersion  spectra  for  0.25  M  (bottom)  and 
2  M  sodium  thiocyanate  with  dispersion  spectra  calculated  for  the  data  of  Fig. 
1  and  a  screening  factor  S  =  0.3  (bottom)  and  0.15  (top). 

Fig.  6.  Same  comparison  as  in  Fig.  5  but  for  0.25  M  (bottom)  and  2  M  sodium 
thiosulfate  and  S  =  0.0  (bottom)  and  -0.1  (top). 
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